Exact periodic and solitonic states in the spinor condensates 
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We propose a method to analytically solve the one-dimensional coupled nonlinear Gross-Pitaevskii 
equations which govern the motion of the spinor Bose-Einstein condensates. In a uniform external 
potential, the Hamiltonian comprises the kinetic energy, the linear and the quadratic Zeeman en- 
ergies. Several classes of exact periodic and solitonic solutions, either in real or in complex forms, 
are obtained for both the F = 1 and F = 2 condensates. These solutions are general that contain 
neither approximations nor constraints on the system parameters. 
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I. INTRODUCTION 



The experimental achievement of Stenger et al. in 
trapping sodium atoms by optical means in 19980, trig- 
gered the study of the magnetism in quantum degenerate 
atomic gases. Since the atom spins are not frozen, the 
direction of the spin can change dynamically through col- 
lisions between the atoms 0-0] • In contrast to the scalar 
gases, spinor gases can host a wide variety of complex 
structures at zero temperature, from spin textures, mag- 
netic crystallization, to fractional vortices et al. 0, @| . In 
the ground state, the symmetry is spontaneously broken 
in several different ways, leading to a number of possible 
phases There exists an interplay between superflu- 

idity and magnetism due to the spin-gauge symmetry. A 
ferromagnetic Bose-Einstein condensates (BECs) spon- 
taneously creates a supercurrent as the spin is locally 
rotated [TJ [l2j]. The study of ultracold spinor is of pri- 
mordial importance in deepening our understanding of 
condensed matter related issues. 

The motion of the dilute spinor condensates is gov- 
erned by the coupled Gross-Pitaevskii equations (GPEs). 
There is a large amount of works that numerically 
solve GPEs [13l4l5| . Analytically, some solitonic solutions 
are obtained by means of variable functions or similar 
transformation for time or spatial modulated coupling 
constants [l6l - [l9T ] . Various approximations are employed 
to study the solitons such as bright and dark solitons 
in the F = 1 spinor BECs (20l - |23| . Exact solutions are 
usually difficult to obtain due to the complexity of the 
coupled nonlinear GPEs. The challenges are two-fold: 
one is the nonlinear density-density interactions, while 
the other is the spin-exchange couplings between the hy- 
perfine states. In our previous publications we have given 
exact solutions for the F = 1 and F = 2 spinor BECs 
for some special cases (2^. [25| . In this paper, we pro- 
pose a general method which simultaneously decouples 
the nonlinear density-density interactions and the spin- 
spin interactions in the GPEs. Classes of the exact so- 



lutions, either in real or in complex forms, are system- 
atically constructed for the Hamiltonian containing the 
linear and quadratic Zeeman energies. The solutions are 
expressed by combinations of the Jacobi elliptical func- 
tions for periodic states or the hyperbolic functions for 
solitonic states. The latter are identified as vector soli- 
tons or scalar solitons, respectively. 

The paper is organized as follows: In Sec. II we de- 
scribed the method and systematically present solutions 
for the spin-1 condensates. In Sec. Ill we give a solution 
to the spin-2 condensates as an example. Section IV con- 
tains a brief summary. 



II. SPIN-1 CONDENSATES 

We are concerned with the quasi-one-dimensional (ID) 
spinor system in a uniform external potential (V(r) = 0). 
In this section, we deal with the F = 1 condensates in 
which the meanfield order parameters are described by 
a macroscopic wavefunction with three hyperfine states 
^ = (ip +1 ,ip , ip-i) T . The Hamiltonian that contains the 
linear and the quadratic Zeeman effects reads 7, 8| 



H 



m— — 1 

Co 



2M 



V 2 +V(r) - pm + qm 2 }ip r , 



a) 



where the spin-polarization vector F = ijj* n F l mn ip n with 

F (i = x, y, z) the spin matrices. The terms with 
coefficients cq and ci describe respectively the spin- 
independent and the spin-dependent binary elastic col- 
lisions in the combined symmetric channels of total spin 
and 2. They are expressed in terms of the s-wave scat- 
tering lengths a and a 2 as: c = ATrh 2 (a + 2a2)/3M 
and C2 = 47r?i 2 (a2 — ao)/3M. p and q are linear 
and quadratic Zeeman coupling coefficients, respectively. 



n-tot = l^i I 2 



V>o| 2 



\ip—i | an d ^( r ) i s the external 



'Corresponding author: yangshijie@tsinghua.org.cn 



potential. 

The dynamical motion of the F = 1 spinor wavefunc- 
tions are governed by idtipm = SH/Sip* n , which are ex- 
plicitly written as the coupled GPEs, 



2 



quasi-lD system. 

Below we choose T% — M — 1 as the units for conve- 
idtipm = [-TTTjdl -pm + qm 2 + c n tot }ip rn (2) nience. By substituting the wavefunction ^(x, t) with 



2M 
l 

c 2 F-FVn,(m= 1,0,-1) 

n=-l 



V>i(M) \ / V^i(a;)e _2(M+Atl)t 
V>oOM) -> Vo(^)e-^ j • (3) 



where cq — c$/2a\ and c 2 = C2/2a\ are the reduced 

coupling constants with a± the transverse width of the we obtain the stationary GPEs as, 



(p + /xOVi = \--dl + (co + c 2 )(|^i| 2 + |Vo| 2 ) + (co - c 2 )|V-i| 2 - p + q]ih + ca^oCi, (4) 
M^o = h^ 2 + (c + c 2 )(|^ 1 | 2 + | 2 ) + col^ol 2 ]^ + 203^5^x^-1, 
(/i - = [-\d 2 x + (c + c 2 )(|^_ 1 | 2 + |^o| 2 ) + (c - c 2 )|Vi| 2 +p + q]iP^ + c 2 ^\. 

I 



Since the chemical potential of each hyperfine state is dif- 
ferent, Eqs.Q are not really stationary equations given 
fii 0. It comprises a "Lamor precession" between the 
hyperfine states. However, the density profile of each 
state are time-invariant so we simply call them the sta- 
tionary states. The periodic boundary conditions, 

1>m(l) = iMO), = C(0). (5) 

is adopted in our calculations. We consider several types 
of real and complex solutions for F = 1 condensates. 

In order to seek the analytical solutions we decou- 
ple spin-spin interactions in the Eqs.Q by requiring 
ip*L m — iipm an d 4>o — ipo- It is directly to check 
that tp-m = V'm corresponds to (partially) spin-polarized 
states (|F| ^ 0) whereas tpl m = — ipm corresponds to 
spin-unpolarized or polar states (|F| =0). On the other 
hand, the nonlinear density-density couplings between 
the hyperfine states are decoupled by making use of the 
properties of the Jacobi elliptical functions or the hy- 
perbolic functions. There are real and complex forms of 
solutions which are explicitly described as follows. 



A. Real solutions 

We first consider the following sn-cn-sn form of solution 
to the nonlinear Eq.Q, 

"01 ( x ) \ I Asn(kx,m) \ 
ipalx) = Dcn(kx,m) , (6) 
ijj-\{x) J y — Asn(kx, m) J 

where sn and cn are the Jacobi elliptical functions of 
modulus m. The period is k — AjK(rn) with j the num- 
ber of periods which will always be set to j = 2 in the 



figures). A and D are the real constants. From the rela- 
tion sn 2 + cn 2 = 1 , one has 

|V-i| 2 = |^i| 2 , |^o| 2 = i5 2 -^|^| 2 . (7) 

By substituting Q into the equations ((4]), we obtain 
three decoupled equations 

AmV'm = ~\^m + TmlVVnlVm, (m. = 0, ±1) (8) 

where the effective chemical potentials /2 m and interact- 
ing constants 7 m are defined as 

Mil = M - 9~ c oD 2 , 
fi Q =H-2coA 2 , 

7 ±1 =2c - Co ^, ^ 
7o = c - 2c ^. 

In order to obtain self-consistent solutions, one should 
take p — — jLii. Therefore, the linear Zeeman energy plays 
the role of balancing the chemical potentials between the 
hyperfine states ipi and ip—i. The effective chemical po- 
tentials (x m and amplitudes A, D can be obtained as, 

fei = ifc 2 (l + m 2 ), 
/i = ifc 2 (l-2m 2 ), 

A 2 = ra 2 k 2 (10) 

7±i ' 
jyi _ m 2 k 2 
70 

From relations © and pH|) we conclude that the effective 
intra-species interactions for hyperfine states ip± should 
be repulsive (7-1-1 > 0) while for ip be attractive (70 < 
0). These impose constraint relations on the values of 
parameters cq, q and [i. 
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FIG. 1: (a) The density profiles of the solution |6]l. Co = 30, 
g = -13.7707, m = 0.4 and ^ = 300. (b) The density profiles 
for the single soliton solution Co = 1.3, k = 2, q = —2 

and fj, — 5. 



FIG. 2: (a) The density profiles of the solution (1121) with Co = 
-10, q = 22.734, m = 0.5 and jit = -200. (b) The density 
profiles for the the single soliton solution corresponding to 
periodic solution (|12[1 with co = —1, k — 3, q — 4.5 and 
fj.= -1.5. 



Figure 1(a) illustrates the density profiles of each hy- 
perfine states for solution (|6|). The parameters are cho- 
sen as c = 30, q = -13.7707, A = 2.0061, D = 2.6704, 
to = 0.4 and [i = 300. This state has vanishing spin- 
polarization |F| = 0. The modulus of the Jacobi elliptical 
functions is a free parameter. As m — ¥ 1, we naturally 
obtain the periodic soliton train solution. The single soli- 
ton solution can be obtained by directly substituting the 
Jacobi elliptical functions with the hyperbolic functions 
in ([6]). Namely, 



V>iO) 
ip (x) 



(Atanh(kx) 
Dsech(kx) 
— Ata,nh(kx) 



(11) 



This solution has been addressed in literatures [lj, Il5j . 
The density profile is displayed in Fig. 1(b), with the 
parameters Co = 1.3 ,k = 2, q — —2, A = 1.6408, 
D = 1.5191 and fj, = 5. It is a dark-bright-dark com- 
posite soliton. 

Other forms of real solutions can also be constructed 
by the same way. For example, we seek the cn-sn-cn form 
of solution to Eq. flD , 

V'lC 2 ') \ / Acn(kx, m) 
i>o(x) = Dsn(kx,m) \ . (12) 
ip-i(x) I \ — Acn(kx, m) 



One has 



l^l 2 , \H 2 = d 2 -^\H 2 - (is) 



Eq.Q are again decoupled into ©. It follows that the ef- 
fective chemical potentials jl m and interacting constants 



A±i =n-q- coD 2 , 
Mo = M — 2c A 2 , 
7±i = 2c - c 

7o 



(14) 



c - 2c -f 



The effective chemical potentials jl m and amplitudes A, 
D are self-consistently calculated as, 



Mil = ifc 2 (l-2TO 2 ), 
fa = i/c 2 (l + TO 2 ), 

_ m 2 fc 2 
_ 7±1 ' 

2^)2 m 2 fc 2 

70 



(15) 



Relations (fT^)) and ([T5|) reveal that the effective intra- 
species interactions for the hyperfine states i/j± should be 
attractive (j±i < 0) while for %p be repulsive (70 > 0). 

Figure 2(a) plots the density profiles of each hyperfine 
states for solution (Tl2)) . The parameters arc Co = —10, 
q = 22.734, A = 3.9602, D = 5.1788, m = 0.5 and 
/i = —200. In this state, the spin is unpolarized |F| = 0. 
Similarly, the single solitonic solution is constructed by 
substitution sn — > tanh and cn — > sech. Figure 2(b) is 
the density profiles of a typical bright-dark-bright soliton 
with Co = -1, k = 3, q = 4.5, A = 2.2913, D = 1.2247 
and )i = —1.5. 
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B. Complex solutions 

We seek complex forms of periodic solution to the non- 
linear Eq.((3]) as, 




/ f(x)e id ^ 
Dsn(kx,m) 



(16) 



where f{x) — y/A + Bcn 2 (kx,m). A, B and D are real 
constants. One has 



\i>-i\ 2 = m\ 



\^ = D 2 (l + ^)-^-\^. (17) 



By substituting (fTTJ) into the coupled GPEs (gj), we 
again obtain the decoupled equations ([5]) with the effec- 
tive chemical potentials and intr-species interaction con- 
stants, 



B 

B) 



Ail = H - q - (co + 2c 2 )^(A + B), 
fr = (J, - 2(c + 2c 2 )(A 
7±i = 2c -(c + 2c 2 )^ 
7o = c - 2(c + 2c 2 ) Jr- 



(18) 



In order to obtain the self-consistent solution, we set 
fix = —p which yields to 



B = 
A = 
D 2 



m k 



2jj ± i-(l-2m 2 )k 2 
37±l 



(19) 



We note that the effective interactions 70 > 0. The 
phase is 



6(x) = 



P(0 



(20) 



where a = ±(2p,±iA 2 - 2j±iA 3 + k 2 AB{\ - to 2 )) i is an 
integral constant. The periodic boundary conditions ([5]) 
require that the amplitude and phase satisfy, respectively, 



/(1) = /(Q), 0(1) - 6(0) = 2j7r x n, 



(21) 



where n is an integer. The periodic condition for the 
phase can be fulfilled by properly adjusting the mod- 
ulus to of the Jacobi elliptical functions. Figure 3(a) 
and (b) display the phase and density profiles of the 
complex solution (|16j) . respectively. The parameters are 
taken as n = 2, c = 34, c 2 = 43, q = 17, A = 1.524, 
B = -0.5795, D = 1.0766, to = 0.82 and fi = 448.6349. 
As to — > 1, it results into the soliton train state. 

The single soliton solution is obtained by substituting 
the Jacobi elliptical functions with the hyperbolic func- 
tions as 




D tanh(fcx) 
f(x)e~ ie ^ 



(22) 
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FIG. 3: The phase profile (a) and the density profile (b) for 
the complex solutions (|16p . The parameters are n — 2, m = 
0.82, g = 17, and p = 448.6349. The phase profile (c) and 
the density profile (d) for the single soliton solution (|22|) with 
n = 6, k = 8.7, g = 5, and /i = 130.69. 



where f(x) = y A + Bsech 2 (kx) . Figure 3(c) and (d) 
plots a typical grey-dark-grey soliton for the solution (|22|) 
with n = 6, Co = 3, c 2 = 4, g = 5, fc = 8.7, A = 
5.2032, B = -2.7032, D = 2.3252 and fx = 130.69. Here 
the phases 6(x) should satisfy the periodic condition ([5]) 
by properly adjusting the width of the soliton fc. Since 
(|F| 7^ 0), this type of composite soliton may be properly 
called the polarized or vector soliton. 

The spin-unpolarized (F = 0) complex solution can be 
constructed as, 




/ f(x)e l9 ^ 

Dcn(kx, to) 
V -f(x)e- ie W 



(23) 



where f(x) = y/A + Bsn 2 (kx, m). A, B and D are real 
constants. One has 



hM 2 = |Vu| 2 , \i> \ 2 = D 2 (l + ^) 



El 

B 



(24) 



The effective chemical potentials and interaction con- 
stants are, 



One obtains 



p,±i=(i-q- c ^-(A + B) 
/j = i u-2c (A + B), 
7±i = 2c - cq%-, 
7o = c - 2c 07 w 



B = 
A = 
D 2 



(25) 



71 



3t±i 
m 2 k 2 
JO ' 



(26) 
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FIG. 4: The phase (a) and the density (b) distributions of 
complex solutions (|23[) for the F = 1 condensates with n — 2, 
m = 0.44, q — —190, if — 4. (c) and (d) are the phase and 
density distributions for the single soliton solutions of (|23[) 
with n = 2, k — 3, q = —8, and /j, — —20. 



FIG. 5: The profiles of phase (a) and the density (b) for F = 2 
BEC solution (I29|l . The parameters are n = 2, m = 0.7, 
q — —53.4153, n — 210. The profiles of phase (c) and density 
(d) for the single soliton solution corresponding to (|29[) with 
n = 5, k = 4.05, q = -8.2012, and n = 200. 



which require that the effective intra-species interactions 
be attractive (70 < 0). 

Figure 4(a) and (b) display the phase and density pro- 
files of the solution (f2U)l . The parameters are n = 2, Co = 
-5, q = -190, A = 4.5347, B = 0.4455, D = 2.7731, 
m = 0.44 and \i = 4. Figure 4(c) and (d) are the corre- 
sponding single soliton solution with parameters n = 2, 
c = -5, q = -8, jfe = 3, A = 0.6781, B = 0.8719, 
D = 1.8825 and \i = —20. It is a typical grey-bright-grey 
composite soliton. In comparison to the polarized soliton 
(|22[) . we may call this spin-unpolarized soliton the polar 
or scalar solitons. 



III. SPIN-2 CONDENSATES 

Our method is also applicable to the F = 2 conden- 
sates. For illustration, we give an example in this section. 
By substituting the wavefunction withjTTl. [25[ 



/ V-2 0M) \ 

W-2(s,t) / 



^i(a;)e- i O+' 11 )* 



(27) 



we obtain the generalized stationary GPEs as 



+ciF T V±i 
„ 1 



CQiitot ± 2ciF z =F 2p- 

C2 



4<# : 



±2 



±1 



c n t ot ± c%F z TP + q}ip±i 



+ Cl (^F T ^ Q + F ± iP ±2 ) - ^=A^ U 

1 \/6 
/»V>0 = [~2^ + c o"-tot]-0o + — c^-F+V'r + B--0-1) 



(28) 



We seek the following form of solution which satisfies 

= and ^=^0, 



/V>2 
!/'() 



\ 



/ 



Csn(fca;, m) 
Dcn{kx, m) 
— Csn(fcx, m) 



/ 



(29) 



where f(x) = \J A + Bsn 2 (kx,m). A, B, C and D are 
real constants. This form of solution has vanishing spin- 
polarization |F| = 0. The equations (|2"9")l can be decou- 
pled and self-consistently solved in the same way as for 
the F — 1 condensates. The linear Zeeman energy satis- 
fies \x\ = —p and fi 2 — 2p so as to balancing the chemical 
potential differences between the hyperfine states. We 
skip the calculation details and just illustrate the results. 

Figure 5(a) and (b) display the phase and density pro- 
files of the solution (|29l) for F — 2 condensates. The 
parameters are n = 2, cq = 10, c 2 = 25, q = —53.4153, 
m = 0.7, A = 1.8, B = 0.6, C = 2.1278, D = 1.7699 and 

= 210. Figure 5(c) and (d) are results for the single 
soliton solution which corresponding to the periodic so- 
lution with n = 5, c = 10, c 2 = 20, q = -8.2012, 
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k = 4.05, A = 3.9, B = 1.3, C = 1.4952, D = 2.429 and 
H = 200. It exhibits a typical grey-dark-bright-dark-grey 
composite soliton structure. 

IV. SUMMARY 

We have systematically solved the one-dimensional 
coupled nonlinear GPEs which govern the motion of the 
spinor BECs exposed in a uniform magnetic field. Both 
periodic and solitonic stationary solutions for the F = 1 
and F — 2 condensates arc constructed. Other forms 



of solutions with different combinations of the Jacobi 
elliptical functions or hyperbolic functions can also be 
obtained in the same way. Our method is general and 
exact, without any approximations or special constraints 
on the system parameters. It may be extended to other 
nonlinear systems such as the coupled nonlinear Klein- 
Gordon equations or the dynamical coupled nonlinear 
Schrodinger equations. 

This work is supported by the funds from the Min- 
istry of Science and Technology of China under Grant 
No.' 2012CB821403. 
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